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ABSTRACT 

A Boltz body of revolution (fineness ratio 7.5:1) was tested in the 
SUMSBS . The effects of sting interference on the drag coefficient of the 
model at zero angle of attack were noted as well as the effects on drag 
coefficient values of boundary layer trips. The drag coefficient values 
were compared with other sources and seemed to show agreement. 

The pressure distribution over the rear of the model with no sting 
interference was investigated including the use of boundary layer trips. 
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1 . I NTR^DUCT ION 

The latest development of the Southampton University Magnetic Suspension 
and Balance System was the completion of a new low speed, octagonal section, 
wind tunnel in 1986. 

A NASA Technical Note, ref. 1, gave drag coefficient values for a Boltz 
body of revolution which appeared to be very low. It was, therefore, decided 
to test a similar body in the SUMSBS low speed wind tunnel. With the body 
suspended magnetically drag data for the true model shape could be obtained, 
but by bringing up a dummy sting to the rear of the model it was hoped to 
simulate sting interference and check the seemingly low drag coefficients. 

An investigation was also carried out into the pressure distribution 
over the rear of the model . As there was no supporting sting to interfere 
with the flow around the base of the model the resulting pressure distribution 
should be more accurate than would be obtained on a sting supported model. 



2 . 


INITIAL DESIGN 


Several areas of design had to be considered for this project. The models 
and associated stings, the optical system, changes to the test section, and 
data analysis programs. 

2 . 1 Models 

The model used in this investigation was the one made for the 
honours project conducted by R. Knight, ref. 2. A set of body surface 
co-ordinates is given for the model in Table 1, and it is shown in 
Appendix 7. The model is comprised of two pieces. The major part 
being the hollow forebody extending to 79% chord, i.e. 79% of the 
total length, the complete model being 184.2mm long. The rear body 
acts as a plug allowing access to the magnets in the model. The 
maximum cross-sectional area, the reference area, is 0.00047 metres 
square . 

Three rear bodies were made, one to accept the dummy sting, one 
to allow pressure plotting over the rear of the model, and one "true" 
shaped to allow testing of a sting free body. 

2.1.1 The rear body designed for the dummy sting. Figure 3, was designed 
to have a cavity in it large enough for the model to be suspended without 
touching the dummy sting. NASA had tested a similarly shape sting supported 
model, ref. 1. Their model was 96 inches in length on a 4 inch diameter 
sting. Scaling this down to the length of the model in this investigation 
of 7.67mm. The cavity in the rear body was bored out to 8mm diameter. 

If testing showed the clearance to be too small the base would be increased 
until satisfactory. The design called for the hollow to be 6mm deep with 
an end face perpendicular to the sides, but a workshop error resulted in 
a M V M shaped end face. 

2.1.2 The rear body designed to allow pressure readings to be taken over 
the rear of the model is shown in Figure 4 and in Appendix 7. The rear 
body was made like a true rear body, but a radial slot was machined out 
as shown. The slot was cut deep enough and wide enough to allow a small 
diameter hollow tube to run along the length of the rear body but beneath 
the surface. The slot was filled in with Araldite, holding the tube in 
position. The Araldite was then smoothed down flush with the surfaces 

of the model. 
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2.1.3 The tube was chosen as 20 cjaucje as this was thought to be the 
smallest diameter that would be easy to bend without causing collapse, 
and reasonably easy to drill into. The tube would emerge on the centre 
line as shown, allowing a rubber tube to be connected to it, connecting 
the model to a manometer. 

2.1.4 Four holes were drilled into the tube at the positions shown in 
Figure 4. The holes were 0.3mm in diameter. Three of the holes would 
be sealed with plasticine at any one time so that pressures could be 
measured at one chordwise station. The plasticine was also used to 
smooth in any indentations in the Araldite . 

2.1.5 The models were made from Aluminium as it is light, easy to 
machine and is magnetically stable. The magnetic core was made up 
of samarium-cobalt magnets as they are magnetically stable and would 
not require continual re-magnetisation as with some other materials. 

Care had to be taken when handling them as they are very brittle. 

2.2 Sti ngs 

Two stings were necessary, one to represent a support sting, and 
the other to act as a pressure take off line. The major design considera 
tions being enough stiffness to reduce vibration at the free end of the 
sting, as little blockage of the tunnel as possible, ease of use, and 
ease of fabrication . 

2.2.1 The stings were designed to be supported on a strut which is shown 
in Figure 5 and is visible in Appendix 7. Rather than have the support 
strut fixed to the inside of the test section it was designed to fix 
to the outside of the test section causing less disruption to the air- 
flow. The strut would then pass through a slot cut into the test section 
With the support strut in place the test section would be too large to 
fit through the axial electromagnets so the strut would have to be fixed 
in place once the test section was positioned in the MSBS. The only part 
of the test section that would take the support strut without having to 
partially dismantle the MSBS was the rear of the test section that 
protrudes from the MSBS. This dictated the length of the sting as it 
would have to reach forwards to where the model was suspended. The axial 
length of the support strut was chosen as 80mm and the width as 10mm in 
order to provide a stable support for the sting. The height of the strut 
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and the shape of the channel to accept the sting were chosen, once the 
sting diameter had been decided on, so that the sting could be positioned 
on the centre line of the test section, 

2.2.2 To provide good stiffness the two stings were designed to fit into 
a common support sleeve which would be attached to the strut. The sleeve 
would reach forwards acting as an outer casing, the required sting could 
then be inserted into the end of the support sleeve and clamped in position. 
The inner diameter of the sleeve was determined by the outer diameter of 
the Support* sting at 5 /16 inch. The outer diameter of the sleeve was 
chosen as 13mm to give a reasonable wall thickness. The sleeve design 

is shown in Figure 7 and in Appendix 7. 

2.2.3 The design of the ’support* sting is shown in Figure 7 and Appendix 

7. The required diameter of the sting was 7.67mm. Brass rod was available 
5 

of diameter /16 inch. This set the inner diameter of the support sleeve. 

As Figure 7 shows the 'support* sting was designed with a separate end 
piece. This had three pressure tappings drilled into it, each connected 
to a 20 gauge tube. This would allow model base pressures to be measured. 
Once the tubes had been soldered to the end piece the end piece was fixed 
to the sting. The first 65mm of the sting was then turned down to 0.3 inch 
or 7.62mm. The tubes were long enough to protrude from the open end of 
the sting, allowing rubber pipes to be connected, and the pressures taken 
off to a manometer. The length of the sting was set at 220mm. This was 
to allow a large range of movement of the front of the sting for changing 
the axial position of the model/sting combination. 

2.2.4 The design of the pressure take off sting is shown in Figure 8 
and Appendix 7. It was designed with a short length of 20 gauge tubing, 
to take the rubber pipe connected to the model, soldered into a larger 
tube to provide stiffness. A short length of tubing was soldered into 
the other end to enable a rubber pipe to take the pressure off to a 
manometer. The larger tube ran through, and was fixed to, a brass rod 
of 5 /16 inch diameter. The whole unit would then fit into the support 
sleeve and be clamped in place. 

2 . 3 O ptics 

The standard set up for the axial position sensing system was to 
have a laser beam fanned out axially, illuminating a long thin detector. 
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The rear of the suspended model would then be in part of the laser 
beam casting a shadow on a portion of the detector. This system would 
not work in this investigation as the Support* sting or pressure take 
off pipe would be in the way of the fanned out laser. 

The two major alternatives were to have a similar system but using 
the front of the model, or a system with the laser fanning out across 
the test section. 

The former would pose problems when trying to launch and retrieve 
the model. It was, therefore, decided to concentrate on the latter 
method. 

2.3.1 There was only a small window in the bottom of the test section, 

so this was taken out and replaced by a perspex window of the same length, 
but greater width, enough to span the flat inner surface of the bottom 
of the test section. Its dimensions can be found in Figure 9. 

2.3.2 It was wanted to illuminate the entire width of the bottom surface 
of the test section with the laser beam. The beam fanning glass rod 
could have been situated in several places but a position just before 
the first mirror was chosen as this would provide the most rigid mounting 
point and would be easier to align with the laser. The laser beam used, 
a lmW helium-neon laser, provides a beam with a Gaussian intensity 
distribution. When the beam is fanned out the beam is less intense at 
its edges than near its centre. Also some intensity will be lost at the 
edges due to greater reflection losses from the two perspex windows. To 
increase the intensity at the edges and make the beam more uniform 

the beam was over expanded and only the central portion directed down 
through the test section. By looking at the action of several diameters 
of glass rod a glass rod with a diameter of 3mm was chosen to fan out 
the laser beam. 

2.3.3 The laser beam would travel from the laser through the glass rod 
and then on to a mirror changing the beam direction so that it would run 
between the top of the test section and the upper electromagnets. It 
would then reach another mirror which would change the beam direction 

so that it passed through the test section onto the detector. 
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Two front reflecting mirrors, 25mm x 16mm, were available and were 
chosen. They were mounted on plastic card and the card was then stuck 
onto the mirror mounts as shown in Figure 10. The card allowed the mirrors 
to be removed from the mounts without damaging the mirrors. 

The mirror mounts were designed to fit on the existing rod and rail 
system within the MSBS . The rods can be moved back and forth above the 
test section, allowing the laser beam to be reflected down through a 
different section of the test section. 

2.3.4 The beam fanning glass rod was mounted on the rod mount shown in 
Figure 10. The mount was then attached to the existing rails in the 
MSBS. 

2.3.5 Initially the existing detector was used, the only change made being 
to mount it transversely. See Section 3 for changes that eventually had 

to be made . 

2 . 4 Manometers 

When running the tunnel two pressures must always be known. Atmospheric 
pressure and the reference dynamic pressure in the tunnel. Atmospheric 
pressure was measured on a standard mercury barometer. In tests with the 
dummy sting the three base pressures must also be measured. This means 
four pressures have to be measured. The tunnel was known to operate with 
a maximum dynamic pressure of 3kPa. Using methylated spirit with a specific 
gravity of 0.83 3kPa is represented by a head of 0.37 metres. The resulting 
manometer height chosen was 0.8 metres which allowed the manometers to be 
used at higher pressures should this be possible in the future. The 
manometers were made by the Chemistry Department's glass blowing service, 
with a bore of 6mm. 

2.4.1 Five such tubes were mounted on a board with a metric graph paper 

scale between the tubes and the board. The board was then mounted on the 

inlet to the tunnel. As the tubes are fairly fragile the board was arranged 
so that it could be removed easily for safer storage. The fifth tube 

al Lowed for any future need to measure five pressures. 
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2 . 5 Data Analysis Programs 

2.5.1 Data Storage 

When running the MSBS the user can store current and position data 

for future analysis. The MSBS runs at 400 Hz, the amount of information 

storeable is set by an array in 4MSBS3 FOR in this case 750 cycles, or 
7 

1 /8 of a second worth of data. Looking at the currents data over a 
second it could be seen that the currents varied considerably. Just 
taking a small portion of data would not reliably represent the currents 
in the electromagnets. Data was, therefore, taken over one second at 
set intervals, 50 cycles total being read. This was arranged by using 
the data file 4MSBS DAT see Appendix 6. When the MSBS was running typing 
lg caused 50 cycles to be stored, the data being taken over a one second 
period. Typing 2g could cause 100 cycles to be stored, the first fifty 
then a seconds break and then the next fifty. 

2.5.2 On stopping the controlling program the data taken could be saved 
into a file of the user's choice for later analysis. 

2.5.3 Program 1 

An existing Fortran IV program was adapted, the new version 4 PROG 1 
FOR being in Appendix 3. On running it would ask for the file name of the 
file in which the current data was stored, along with the room temperature 
and pressure in Celsius and kPa respectively. It then asks the user how 
many blocks of data are to be read from the data file and whether the 
output is to the screen ^ or the printer When running the MSBS 

typing lg results in 50 cycles worth of data being stored. The data is 
read over a second as five blocks of ten cycles. Hence to read in lg 
worth of current data the user must ask for five blocks to be read in, 
ten blocks for 2g, etc. 

2.5.4 The program reads in the required number of blocks and calculates 
the average of the sum of the axial currents of each cycle, and also the 
standard deviation from that average. Using Chauvenet's criterion for 
error rejection any 'bad' cycles are rejected the new average currents 
and positions being calculated. 

2.5.5 The program then asks for the date when the data was taken. The 
date to be inputted in a coded form. If the wind was off and the data 


7 


is the first set then the format is "START date". If the wind was off 
and the data is the last set then the format is "STOP date". If one of 
these is used the progrm asks if the user wishes to continue. If the 
data is not the first or last set then a different format must be used. 

If the dummy sting was present when the data was taken the format is 
"YS date", the program then goes on to ask for the three base pressures 
and the reference dynamic pressure in kPa, and then on to ask if the user 
wishes to continue. If the dummy sting was not present the format is "NO 
date", the program then asks only for the reference dynamic pressure, and 
then on to ask if the user wishes to carry on and analyse the next block (s). 

2.5.6 Once all the data has been read in the averaged current data for 
each block along with average position data, pressures and temperature, 
and the date are saved in the file 4PR0G1 STO which is used in the next 
program. The data can also be added to a master file 4CURR DAT with the 
averaged currents of all the runs to date stored in it. 

2.5.7 The data is saved with the "START date" and "STOP date" being 
saved first. This means that when the file is read in the next program 
the current offsets can be changed before the wind on currents are read 
in. 

2.5.8 Program 2 

A listing of 4PROG2 FOR can be found in Appendix 4. After running 
4 PROG 1 FOR this program is run. It reads in the contents of 4PR0G1 STO 
one line at a time. If the data is a "START date" or "STOP date" then 
the respective current offset is reset and then next line is read in. 

If the date is not a "START date" or "STOP date" the program goes on to 
convert the axial currents to an axial force using a factor obtained from 
the current force calibration of section 4.2. The test section dynamic 
pressure is then calculated by multiplying the reference dynamic pressure 
by a value obtained from the pressure calibration in section 4.1. The 
density of the air and its kinematic viscosity are then calculated. It 
was assumed that the relative kinematic viscosity follows the relative 
absolute temperature to the power 0.8. 

2.5.9 If the date signifies that the sting was not present the base 
force is set to zero, otherwise the baseforce is calculated. Knowing the 
three measured base pressures and the radii at which they were measured 
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the program assumes that the pressure varies as a quadratic function of 
radius. Calculating this function , see Appendix 2, the pressures acting 
on the base area are integrated up to give the acting baseforce. 

2.5.10 The base force and the force provided by the axial magnets act 
in the same direction. The sum of the two is equal to the drag force 

on the model. Knowing the drag force, the dynamic pressure, and the 

-3 +2 

reference area, 0.47.10 m the maximum cross sectional area, the drag 
coefficient is calculated. The test section free stream velocity, Mach 
number, and Reynolds* number, based on the length of the model, are then 
calculated. 

2.5.11 The drag coefficient, Reynolds* number, free stream velocity, free 
stream Mach number, free stream dynamic pressure, total axial force (equal 
to the total drag force acting on the model) , the baseforce, and the date 
are then written to the output file 4PROG2 STO. 

2.5.12 The program then returns to read in the next line of data from 

4 PROG 1 STO. If the line read in is the *'END FILE" line the program asks 
if the user wants the data (Cd etc.) to be written to the master file 
4CDRAG DAT. The program then saves the necessary files and then stops. 



3. 


PRELIMINARY TESTING AND REDESIGN 


3-1 As stated in Section 2.3.5 the axial position detecor was initially 
mounted transversely. A red ’Cokin' f i£ter was used to cut out as much 
ambient light as possible and also to protect the detector. 

3.1.1 The model was suspended in the MSBS. It was found that the model 
oscillated axially about ± 5mm from an initial position. With the particular 
shape of the rear of the body being tested, i.e. a low curvature until 

near the end of the model, an axial movement only produces a small change 
in the shadow size and hence change in the output of the detector. With 
the model out of the MSBS the output of the detector fluctuated, due to 
variations in the intensity of the laser beam with time. This could be 
regarded as noise, and any change in output due to the model moving as 
the signal wanted. The axial oscillation would be reduced by increasing 
the signal to noise ratio. 

3.1.2 This was initially performed by blanking off parts of the detector 
which previously would have always been illuminated regardless of the 
model's axial position. The oscillation was greatly reduced, but at times 
when the laser was fluctuating badly was still too high. 

3.1.3 Several options were open; use a more consistent laser, postion 
the laser at the front of the model with an axially fanned laser beam, 
use a digital sensor which would not interpret the slight changes in 
intensity as a position change, or use another detector to provide a 
reference signal indicating any changes in the intensity of the laser 
beam. 

3.1.4 The first option was not possible. The second was discarded for 
the reasons given in Section 2.3. The third would pose a problem as the 
existing digital sensors are only about one inch in length and would 
have either required a collimated laser beam, or would have only been 
able to pick up the laser/shadow boundary on one side of the model if 
the beam was diverging as it had been set up to do. Only picking up the 
one boundary would pose problems as any yaw of the model would be inter- 
prted as a large axial position change, if two boundaries were used this 
effect would nearly cancel out. The fourth option of using a reference 
sensor appeared to be the most likely solution given the time available. 
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3.1.5 A small, 6mm x 6mm photodetector was fixed in the same housing 

as the main detector. The reference detector was connected to an 

a 

amplifier and anlogue to digital converter using the "roll" channel 
which would have been unused in this investigation. The voltage from 
the reference sensor was present on channel one, and the voltage from 
the main detector on channel zero. 

3.1.6 An existing Fortran IV program was modified to repeatedly read 

in both voltages and display them along with the changes from the previous 
values. The changes did not keep in step. The program was altered to 
display the changes from the first values read in. This gave much better 
results with the changes keeping in step for large fluctuations from the 
original values. For smaller changes of the reference detector the 
results were not so good. A threshold was built in so that only changes 
greater than a set value would result in a change being made to the value 
obtained from the main detector. The next step was to incorporate a 
MACRO version of theprogram into the control program, 4MSBS3 MAC, used 
to control the MSBS. 

3.1.7 The final version of 4MSBS3 MAC is in Appendix 5. When the program 
is run the threshold value and scaling value are set at zero. These 
values could be altered while the program is running by typing in the 

new value along with a letter signifying the parameter to be changed, 
and then typing "l?" which initiates the changes. The reference outputs 
of the main sensor (setting the model's axial position) and of the 
reference sensor can be changed by the same method. The scaling factor 
used in the program has to be inputted as an integer. To allow a factor 
between one and two (which appeared to give the best results) the required 
factor is inputted in thousands i.e. 1500 for 1.5. For example to set a 
factor of 1.75, a threshold of 3, to move the model to a position such 
that the main sensor output is 2550 and to reset the reference sensor 
reference value to the instantaneous ready type 1750F 3 t 2550M IS 1?. 
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4 . CALIBRATIONS 

Two calibrations were necessary. One to relate the reference dynamic 
pressure to the dynamic pressure in the test section, and the other to relate 
the applied axial current to the acting drag force. 

4 • * Pressure Calibration 

The reference dynamic pressure was measured by a pitot tube and a 
static pressure tapping just aft of the contraction. The total pressure 
of the air flow was assumed constant along the tunnel. A static pressure 
tapping already existed in the test section. Two manometers were used 
to show the true dynamic pressure and the reference dynamic pressure. 

Also the difference between the two static pressures could be calculated. 

A plot of the two dynamic pressures is given in Figure 11. The relationship 
q^H^E = ^ref‘ 1,017 Can determined f^om the plot. 

4.2 Force-C urren t Calibration 

A static force-current calibration was used. An axial force was 
applied to the model, the resulting currents recorded and a plot of the 
results drawn up, see Figure 12. The force was applied by fixing one 
end of a length of cotton to the near of the model, the other running 
over a pulley to a pan with weights on it. The plot has a gradient of 
0.044 Newtons per amp, where the current is the sum of the axial currents. 


i 
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TESTING 


5. 

The testing was carried out in two blocks. Firstly the drag data was 
taken with and without the dummy sting, with and without boundary layer trips, 
and secondly an investigation into the pressures over the rear of the model 
with and without boundary layer trips was undertaken. 

5. 1 Dr ag Investigation 

The model was suspended and moved into the position it occupied in 
the force-current calibration with the commands 280-m 300y 40z 2550M IS 
2 t 1?. Any roll was taken out by holding the model so that it stopped 
oscillating and by releasing it carefully. The tunnel intake was then 
wheeled up to the test section, aligned using the aligning dowels and 
clamped to the test section. If the model was oscillating axially the 
oscillation was reduced as much as possible by inputting a series of 
scaling factors and then using what appeared to be the best one, e.g. 

400F 1?. 

5.1.1 Current data was recorded with the wind off. The wind was then 
turned on with the tunnel speed control set to either 3iV or 6V (the 
maximum is 10 Volts) to avoid low speed resonances in the fan unit. 

The required speed was then set and the flow allowed to stabilise for 
around a minute. Current data was taken and the necessary pressures 
noted down. Once all the wind on data had been taken the tunnel control 
was set to either 3iv or 6V the tunnel allowed to slow down, and then 
the fan unit switched off. Once the fan had stopped spinning and there 
was no, or negligible, flow through the tunnel a set of wind off data 
was taken. 

5.1.2 The tunnel intake was removed from the front of the test section, 
and the model removed from suspension. The escape character was then 
pressed and the data saved into a file. The electromagnets were then 
switched off. 

5.1.3 It had been planned that the sting should reach into the base 
cavity. In practice it was found that this was easy to do but that the 
model kept touching the sides of the sting due to " sloppiness" in the 
control system. It was, therefore, decided to suspend the model just 
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in front of the sting rather than to increase the diameter of the cavity. 
The change in model position had no noticeable effect on the base pressures 
so it was felt to be valid. Also the change in position meant that the 
integrators could be used to improve the position holding of the model. 

5.1.4 The data was then analysed using the programs 4PR0G1 FOR and 
4PROG2 FOR with the calculated drag coefficients etc appearing in the 
file 4PROG2 STO. 

5.1.5 Data was taken for several conditions. The sting-less model was 
tested with no boundary layer tip, and boundary layer tips at 40%, 60%, 

75% and 85% of the chord of the model. The model with the dummy sting 
was tested without a boundary layer tip, and with a trip at 85% of the 
model’s chord. The trips were made from two strands of copper wire 
twisted together and then soldered in a loop. The resulting trips were 
0.012 inches in diameter. If necessary they were stuck to the model 
with a 1 super glue 1 . 

5.2 Pressure Investigation 

The rear body with the pressure take off ’tube 1 had four static 
holes drilled through the araldite into the pressure tube. For testing 
three holes were smoothed over with plasticene so that only one hole 
was open. 

5.2.1 The pressure take off sting, see Section 2.2.4, had connected 
at its rear end around four feet of small bore rubber tubing. The 
other end of the tubing passed through the support sleeve of the sting, 
out through the top of the test section, and was connected via an adapter 
to a larger bore pipe. The larger bore pipe was connected to one side 

of one of the manometer tubes, the other side being connected to the 
reference static pressure. 

5.2.2 To launch the model required two people. The sting would be 
unclamped from the support sleeve and pulled forwards to the front of 
the test section. It was then connected to the model via a short length 
of tubing. The tubing connected to the other end of the sting was then 
pulled so that the sting and the model moved down the tunnel. The model 
was supported on a makeshift stand inside the test section while the 
sting was clampled into the support sleeve. Care had to be taken to 
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ensure that when the model was moved forwards to the limit of the short 

* 

tubing it remaind illuminated by the laser beam. Holding the model in 

A 

one hand the stand was removed from the test section. Making sure that 
the hand holding the front of the model was not in the way of any light 
beams and that the model was central in the test section the controlling 
program was started by the other person, and the electromagnets switched 
on. If the MSBS did not pick up the model successfully the magnets were 
switched off. They were then switched on again. With practice the model 
launched first time. 

5.2.3 The tunnel intake was then fitted, as in Section 5.1 and the tunnel 
turned on at 3iV or 6V. The required speed was set, the flow allowed to 
stabilise and the reference dynamic pressure and the difference between 
the pressure from the open static tapping on the model and the reference 
static pressure noted. The tunnel speed was then changed and the new 
pressures noted. Once the required data had been taken the tunnel was 

slowed down and then switched off. The tunnel intake was then removed. 

r 

One person would hold the front of the model taking care not to interupt 
any of the light beams. The other person would then switch off the 
electromagnets and stopping the controlling program. The stand was then 
used to support the model . 

5.2.4 To make an alteration to the model, such as to put on a boundary 
layer trip, or open up a different static tapping, the sting was unclamped 
from its support sleeve and pulled forwards along with the model to the 
front of the test section. The model could then be disconnected from the 
sting and any alterations made. 

5.2.5 Each of the four static tappings was used at the six speed settings 
used in the drag investigation with no boundary layer trip, and with 
boundary layer trips at 75% and 85% of the model's chord. The testing 
order was to seat three tappings, test with no trip at all the speed 
settings, then with the trip at 75% chord, and finally with the trip at 
85% chord. The next tapping was then opened up, the previous one being 
sealed. 

5.2.6 The pressure coefficients were calculated using the method in 
Appendix 1, the results being shown in Figures 15, 16 and 17 and Table 2. 
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6 . 


RESULTS AND DISCUSSION 


The investigation was carried out in two parts. The drag investigation 
was carried out first, and then the pressure distribution investigation. The 
results, however, have to be viewed as a whole with the drag results directly 
related to the pressure distributions. 

6 . 1 Comparison of Sting Free and Sting Based Models* Cd Values 

Figure 13 shows the variation of drag coefficient with Reynolds number. 
It can be seen that the Cd values for the model without a sting and with 
no artificial boundary layer trip are higher than the values for the model 
with a sting and without a trip. On the model without a sting the flow 
separates near the base due to the strongly adverse pressure gradient 
around the base, and a wake is formed. With the sting in position the 
rear body-sting shape is not curved so sharply, the pressure gradient 
is still adverse but much less so, and the flow is able to follow the 
shape of the rear body-sting. Separation may occur due to the radial 
distance between the end of the body and the sting, but it would be less 
severe than on the sting free model. As the flow around the base is 
more favourable the pressure drag is lower for the rear body-sting 
combination, and hence the drag coefficient is reduced. 

6.1.1 The general trend of the trip-less models is that the Cd values 
fall as the Reynolds* number increases. This can be explained by the 
relative thinning of the boundary layer with increasing Reynolds* number. 

The tests were carried out over a low Mach number range (Max. 0.21) and 
over a low Reynolds* number range (Max. 8.7.10 5 ), the fall in Cd values 
was therefore expected. The fall is more consistent for the model with 
the sting than without the sting. 

6.1.2 A possible reason for this could be that the separation of the 
boundary layer from the sting free model is unsteady, with nothing to 
fix the separation. Or the variation could be due to errors in the 
testing procedure. The tests on the sting free model were the first 
tests, and the testing procedure was subsequently improved. As described 
in Section 2.5.1 current data was taken in small batches over a second 

in an effort to obtain a better average of the currents. However, this 
system was not used initially. The currents for the sting free model 
with no boundary layer trip being taken continuously over an eighth of 
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a second, and hence more likely to pick up temporary fluctuations in 
the axial currents. This would follow through to the drag data. 

6 . 2 Comparison of Resul t s with other Sources 

Figure 14 shows a comparison of the two trip free models' drag results 
with those obtained in references 1 and 3. It must be born in mind that 
the results obtained in this report are not corrected for wind tunnel 
interference effects, whereas the values from references 1 and 3 are. 
References 5 and 6 were consulted but did not give a useable method for 
a closed regular octagonal test section with a body of large chord relative 
to the test section height. Opposing flats in the test section are seven 
inches apart. The expected correction to the Cd values would be a decrease 
in the region of one to two percent. 

6.2.1 Looking at the two solid line plots shows the results for sting 

based models. The results from ref. 1 are over a Reynolds' number range 

an order higher than the range attainable in the SUMSBS wind tunnel. They 

were taken at a constant Mach number of 0.1. The plot shows that as the 

Reynolds' number increases from 5.10 the Cd values increase. The lower 

Reynolds' number points show the curve to be levelling off. The Cd values 

obtained in the SUMSBS wind tunnel fall with increasing Reynolds' number. 

The two sets look as if they may join up. It would be interesting to take 

6 6 

some data in the range 1.10 < Re < 4 . 10 to see if the sets can be linked. 

The Cd values in ref. 1 are given as based on wetted area. They are shown 
in Figure 14 as being based on maximum frontal area, as are all the other 
Cd values. 

6.2.2 The broken lines show the sting free models, from this investigation 
and from ref. 3 for an ellipsoid of fineness ratio 8:1. Both sets fall with 
Reynolds* number. The ellipsoid has a lower drag coefficient than the 
Boltz body used in this investigation, by around ten percent. Reference 3 
suggests that as the ellipsoid's fineness ratio decreases the Cd value 
should fall, with the optimum ratio being 4;1. The Boltz body used has 

a fineness ratio of 7.5:1, the difference in shape must therefore account 
for the higher drag, and it suggests that for a given fineness ratio the 
Boltz body has over ten percent higher drag. Correcting the Boltz body 
for wind tunnel interference would only be expected to reduce the Cd values 
by around one to two percent. 
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6 . 3 The Effects of Trip ping the Boundary Lay er 


Referring to Figure 15 the effect of using the boundary layer trips 
can clearly be seen. The trip alters the Cd value by several means. There 
is a change in Cd due to the drag of the wire itself , although this should 
be small relative to the total drag of the model. If the boundary layer 
is laminar when it reaches the trip it separates, to reattach further 
downstream as a turbulent boundary layer. The turbulent boundary layer 
produces a higher skin friction than a laminar boundary layer so the Cd 
will rise. However, the turbulent boundary layer is more energetic and 
is more able to withstand an adverse pressure gradient and the separation 
occurs further downstream. This results in a smaller wake and more 
pressure recovery over the base of the model, and the Cd falls. The total 
change in Cd depends on the position of the trip, too far forward and the 
skin friction outweighs the pressure recovery and the Cd value rises. 

6.3.1 With the tip at 40% chord on the sting free model the Cd values 
increase dramatically over the no trip model. Initially the Cd rises to 
a peak and then falls off as Reynolds' number increases. This could be 
due to the pressure recovery over the base only being large enough to 
outweigh the increase in skin friction at higher Reynolds' numbers. 

6.3.2 With the trip at 60% chord the behaviour is as for the trip at 
40%, but not so great. 

6.3.3 With the trip at 75% chord on the sting free model the Cd values 
are lower than the values for the no trip model. The trip is now acting 
beneficially, the pressure recovery outweighing the increase in skin 
friction aft of the trip. The Cd values fall with increasing Reynolds' 
number, but with an extra low Cd at a Reynolds' number of 7.2.10^. 

6.3.4 With the trip at 86% chord on the sting free model the same thing 
happens at a Reynolds' number of 7.2. 10^. The extra low value is not 
such a large variation from the other points as for the 75% trip plot. 

6.3.5 From Figure 15 it can be seen that at the rearmost pressure tapping, 

at 97% chord, the Cp values drop slightly (become less positive) at a 

5 

Reynolds' number of around 7.2.10 . This would imply that the normal 
force at this position was decreasing, i.e. that the drag should increase. 
The Cp values for the tapping at 91% chord do not show a major change at 
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this Reynolds' number. These results are anomalous with the dip in Cd 
values at this Reynolds' number. This suggests that the dip in Cd values 
is due to changes in the pressure distribution over parts of the model 
that were not investigated, i.e. aft of 97% chord. Another possibility 
is discussed in Section 6.4.4. 

6.3.6 The plot for the sting based model with a trip at 86% chord also 
shows strange behaviour at a Reynolds' number of 7.2. 10^, but in this 
case it is an increase in Cd value. There is no Cp information available 
to explain this. 

6.3.7 The plot for the sting free model with a trip at 86% chord shows 
that at low Reynolds' numbers the Cd value is higher than the no trip 
values. As the Reynolds' number increases the Cd values drop below the 
no trip values. This low Reynolds' number behaviour does not occur with 
the trip at 75% chord, which has Cd values slightly lower than the no 
trip values, but it does occur with the sting based model with a trip 

at 86% chord. 

6.3.8 From Figure 16 it can be seen that at the lowest Reynolds' number 

of 3.6.10^ the Cp values with the trip at 75% chord only differ appreciably 
from the no trip values over the last nine percent of the chord, or maybe 
less if a different line is drawn through the points. As the axial position 
increases the Cp values become less negative until at 93% chord they become 
positive values, i.e. the pressures reduce the drag relative to the no trip 
model. The difference between the Cp at the 97% chord tapping is large, 

0.3 with the trip as against -0.5 without the trip. With the trip at 86% 
chord there is a major change in Cp behaviour with axial position. The 
overall trend is still for the Cp values to become less negative with 
axial position, but until 94% chord, or more, the Cp values are more 
negative than the values for the no trip model. It seems therefore that 
the more negative Cp values, 86% chord trip relative to no trip, from 
80% chord to 94% outweighs the less negative values over the last 6% 
of the chord and the Cd is actually increased relative to the no trip 
model . 

e 

6.3.9 With the trip at 75% chord the bounday layer reattachs after 

A 

the trip, but with the trip at 86% and the model curvature much higher 
it may be that the boundary layer does not reattach. 
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6.3.10 As the flow is subsonic pressure information can be transraiteed 
upstream and the separation could affect axial positions upstream of the 
trip. 

6.3.11 The Cp values for the trip at 86% chord cannot be explained by 

a leak in the pipes carrying the pressure to the manometer as the Cp 

5 

values for the Reynolds* number of 8.8.10 do not show any unusual 
behaviour , Figure 17, the two sets of data being taken at the same time, 
see Section 5.2.5. 

6.3.12 Figure 17 shows the pressure variation over the rear of the 

5 

model at a Reynolds* number of 8.8.10 . The plot for the model with no 
trip shows the same trend as at the lower Reynolds' number, but the Cp 
values at 97% chord is less negative and therefore results in less drag. 
The plot for the trip at 86% chord shows the Cp values becoming less 
negative and finally positive reaching a Cp of +0.08 at 97% chord, 
accounting for the decrease in Cd over the no trip cd value. With the 
trip at 75% chord pressure recovery occurs aft of 91% chord, but between 
81% and 91% chord the Cp values are more negative than the no trip values. 

6.3.13 Figure 13 shows that at the higher Reynolds' number the trip at 
86% chord on the sting based model results in lower Cd values than the 
no trip sting based model. Possible reasons could be due to a limited 
pressure recovery over the remainder of the model's rear surface, or 
that the change in flow between the model and the sting is accomplished 
more easily and clearly with less disruption being transmitted upstream. 

6 . 4 Sources of Erro r 

The model used was representative of the Boltz body as defined in 
Table 1. However the model was not well made, deviating from the given 
co-ordinates. The model used in ref. 1 on the other hand was of extremely 
high quality. Also the construction of the model with a main forebody 
and a smaller rear body meant that the two sections had to be joined 
together. The joint between the two was not perfectly flush. It was 
feared that the joint may initiate transition, but this did not appear 
to be the case. The joint was at 19 % chord. Figure 13 shows that with 
a trip at 75% chord there is a large drop in Cd values over the 'no trip* 
values. If transition was occurring at the joint the trip at 75% chord 
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would produce higher Cd values than the 'no trip' values due to the 
extra skin friction and wire drag. Therefore it seems that the joint 
does not initiate transition. 

6.4.1 By looking at the currents in the axial electromagnets with the 
model in suspension and the wind off it could be seen that the axial 
currents varied by around ± 20% from a mean value. This variation also 
incorporated the variation due to the axial position sensing system 
picking up some of the fluctuations in intensity of the laser beam. 

This variation in axial current follows through directly into the Cd 
values. It is expected that the residual variation in axial currents 
would stay the same as the current increased, so the percentage change 
would drop. This appeared to happen with the scatter in Cd values being 
higher at lower Reynolds' number. 

6.4.2 Each point on Figure 13 is the average of fiteen values. Each of 
the fifteen values is made up from fifty data points taken in five batches 
of ten data points. In this way it was hoped that the average value would 
be representative of the true value. With the average being taken of 

the fifteen values if one of those values was increased by twenty percent 
the average would change by less than two percent. 

6.4.3 The pressure investigation was carried out with a small diameter 
(20 gauge) tube protruding from the rear of the model, see Section 2.1.2. 
The rearmost pressure tapping was at 97% chord and it is felt that the 
pressures measured there would only be slightly affected by the tube. 

6.4.4 The flow through the test section has not yet been surveyed. The 
flow seemed steady, but it is interesting to note that the flow was not 
as steady with the fan voltage at 8V which corresponds to a Reynolds' 
number of 7.2. 10 5 . This was the Reynolds' number at which the kinks 
appeared in the drag coefficient vs Reynolds' number plot see Figure 

13 and Section 6.3.3 to 6.3.6 inclusive. The models were positioned 
in the centre of the test section. 


21 


7. 


CONCLUSIONS 


The drag coefficient values obtained suggest that the seemingly low values 
attained in ref. 1 could be correct, and due to the change in base flow due 
to sting interference, and to Reynolds' number effects. 

The general trend of drag coefficient with Reynolds' number, over the 
Reynolds' number range attainable in the SUMSBS low speed wind tunnel, is that 
the drag coefficient values fall with increasing Reynolds' number. 

The drag coefficient values for sting free model are generally 15% higher 
than the values for the sting based model. This shows the importance of allowing 
for sting interference effects, and the benefits of an MSBS. 

The use of boundary layer trips has a large effect on drag coefficient 
values, and hence when testing obtaining a similar transition point to full 
scale is crucial. 

The MSB S allowed pressure plotting over most of the rear of the model. 

The results showed the expected trends, although data from more axial position 
would have been useful. The method worked easily, but there was room for 
improvement. 
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8 . RECOMMENDATIONS 


8.1 The present MSBS and wind tunnel could be improved by: 

1. The further development and use of a digital axial position 
sensing system which would not be affected by the small fluctua- 
tions in the intensity of the lower beam. 

2. Combining the fan speed control and on/off switch in one unit 
next to the wind tunnel . 

3. The control system program be optimised to provide better 
model control. 

4. The elimination of the low speed fan unit resonances. 

5. Fitting a trap door in the diffuser allowing easy retrieval 
of any model flyaways. The model was only lost twice, both when 
the fan was turned straight from full on to off and even then 
the model stayed in the test section. 

6. The fitting of a seal between the test section and the air 
intake. 

8.2 Suggestions for further investigations are: 

1. To carry out an accurate velocity traverse of the test section. 

2. To carry out sting interference tests with the model at various 
angles of attack. 

3. To extend the Cp investigation further to the rear of the model/ 
and also to check the values obtained in this report. 

4. The possible development of a solid-state on-board pressure 
sensor able to transmit data to a receiver without a physical link. 
This would allow testing right up to the trailing edge of the model, 
with no interference. 
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Fig 11 Dynamic Pressure Calibrati on 
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Fig 11 Axial Force-Current Calibration 
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Fig lfe Variation of Cp with Axial Position 
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Table 2 


pleasured values of Cp and Reynold** No. 


tour static tap pi nqs 


Chord 31% 

Cp _Re 
1E5 


No trip:- 

-.061 3.67 
-.064 5.77 
-.065 6.55 
-.053 6.32 
-.063 7.31 
-.062 7.60 
-.053 3.03 
-.053 3.22 



-.043 3.61 
-.033 5.77 
-.046 6.55 
-.045 6.34 
-.042 7.23 
-.045 7.54 
-.044 3.03 


chord i 
67 I -, 
63 I - 


Trip at 75% 
-.061 3. 
-.047 5. 
-.051 6. 
- . 050 6 . 
-.043 7. 
-.048 7. 
-.044 3. 
-.052 8. 


Trip at 36% chord:- 

-.061 3.67 -.064 3.53 
-.072 5.73 -.053 5.75 
-.072 6.51 -.051 6.52 
-.063 6.73 -.053 6.73 
-.070 7.23 -.053 7.30 
—.063 7.62 —.051 7.64 
— .061 y • 0 0 — .051 3.03 
-.064 8.73 -.051 8.76 


31% 

Cp Re 

37% 

Cp Re 


1E5 


1E5 

-.013 

3 . 65 

-.052 

3.48 

-.016 

5.53 

-.041 

5.63 

-.023 

6.44 

-.027 

6.43 

-.023 

6 . 78 

- 

- 

-.022 

7.23 

-.033 

7.35 

-.021 

7.73 

- 

- 

-.025 

3.04 

-.031 

3.03 


3.78 

-.023 

3.76 

i 

-.014 

3.60 

.032 

3.60 

-.016 

5.63 

.050 

5.53 

-.023 

6.44 

.057 

6.33 

-.024 

6.76 

- 

- 

-.022 

7.23 

.055 

7.27 

-.022 

7 . 63 

- 

- 

-.025 

3.02 

.056 

7.35 

-.023 

8.75 

.057 

8.73 

-.048 

3 . 60 

.001 

3.64 

-.022 

5 . 65 

.081 

5 . 66 

-.013 

6.43 

.081 

6.43 

-.013 

6.73 

- 

- 

-.014 

7 . 25 

.073 

7.41 

-.015 

7 . 66 

- 

- 

-.012 

7.33 

.073 

3.13 

-.013 

8.73 

.083 

3.34 
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APPENDIX I 


CALCULATION OF PRESSURE COEFFICIENTS 


The definition of pressure coefficient is 


Cp * 


P -P 
m oo 


oo 


where 


p^ is the measured pressure 

P oq is the free stream static pressure 

q is the free stream dynamic pressure. 

The tunnel is set up so that the static pressure at the model, p , can 

m 

be measured referenced to the reference static pressure, p^, measured just 
after the contraction. 


The measured difference is p -p 

m r 


P ^P 
m oo 


(P "P > - (P -H) + (P -H) 

m r oo r 


where H is the total pressure of the free stream. 


Now, q is related to the reference dynamic pressure, see Section 
oo 

4.1, by the relationship 

q - q . 1.017 = (H-p ) . 1.017 

*oo m REF *r 


Cp 


P “P- 
m r 

(H~P r ) 1 TO 1 7 


+ 1 + 


p -H 
r o 


(H-p r ) 1.017 


. . Cp = 1 + 


p -p 
m r 


1 


(H-p ) 1.017 


1.017 


Hence the pressure coefficient can be calculated from the two known pressure 

differences p -p and H-p . 

m r r 
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APPENDIX II 


CALCULATION OF BASEFORCE 


The baseforce acting on the model with the dummy sting has to be measured. 
By measuring three pressures at known positions a pressure distribution is 
found and used to calculate the baseforce. 

Assuming a quadratic pressure distribution around the centre line of 

2 

the base cavity of the form P(r) = a + br + cr . 




The 

element of area 

is 

S A (r ) 

= 2tt r Sr 



. The 

element of force 

is 

SF(r) 

= P(r) . S A (r ) 

= 2 rr (ar 

w 2 3. . 

+ br + cr ) dr 





R 



.The 

total baseforce 

is 

F = 

\ 2 rr (ar + 

) 

br 2 + cr 3 ) 

dr 




F = 

o 

, f ar 2 

2tt [ — . 

br 3 

3 + 

4,' 

" lo 




F ~ 

2 " i -r ♦ 

bR 3 

3 + 

} 


where R is the maximum radius of the cavity. 

The coefficients a, b and c are found by considering the three known pressures 
at the three known radii. See the listing of 4PROG2 FOR in Appendix IV. 
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APPENDIX III 


4PR0G1 . FOR 


C 

C 

C 

C 

C 

C 

200 

210 

805 

C 

822 

900 

C 

876 

877 

101 

102 

C 
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PROGRAM 4PR0G1.F0R 

A. W. NEWCOMB 

DIMENSION I DAT A < 1 2000 > f NAME < 6 > f DATES ( 7 ) 

DIMENSION AUERS ( 71 f 8 ) f ARRAY < 50' 14) 

DIMENSION IFLAG2<50) fSUM< 14) f AV< 14) fAVBACH<8> 

INTEGER AXAVfHfAfP 
BYTE KEY 

IC0UNT=0 

IGRP=0 

TYPE 200 

FORMAT (//'$ FILENAME OF DATA SOURCE — ') 

READ(5r210 ) (NAME(I) fI=1f6) 

FORMAT (6A2) 

WRITE<5f805) 

FORMAT ( ' RESULTING DATA TO BE STORED IN 4PR0G1.ST0 ') 
WRITE(5f822) 

FORMAT ( 7 INPUT ROOM TEMP CELSIUS AND PRESSURE KPA ') 
READ<5f*)RTCfRP 

OPEN ( UNI T=1 f NAME=NAME f TYPE* ' OLD 7 f FORM* ' UNFORMATTED 7 f 
1REC0RDSIZE*! ) 

OPEN < UNIT*2 f NAME= 7 4PR0G1 * TEM 7 f TYPE* 7 SCRATCH 7 f FORM* 7 FORMATTED 7 ) 

OPEN ( UN I T=3 f NAME= 7 4PR0G1 ♦ STO 7 f TYPE* 7 NEW 7 f FORM* 7 FORMATTED 7 > 

CONTINUE 

AASUM=0 » 0 

AXSUM-0 *0 

FASUM=0 * 0 

AAAU=0 • 0 

AXAV=0 

FAAV=0 *0 


FORMAT ( 7 HOW MANY BLOCKS OF DATA TO BE ANALYSED -- 7 ) 
WRITE < 5 f 876 ) 

READ (5f 877) I SETS 
FORMAT (21 ) 

NUMBER=10 
NR=NUMBER 
ND-NUMBER*16 
TYPE 101 

FORMAT ( / 7 ♦ TERMINAL(5) OR PRINTER(6) 7 ) 

READ ( 5f 102 ) I PLACE 
FORMAT (ID 
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953 CONTINUE 

DO 10 1 = 1 f ND 

10 READ(l) I DATA ( I ) 

C 

XASUM=0* 0 
XFSUM=0 ♦ 0 
WRITEdPLACEf 20) 

20 FORMAT < 30X f 15HCURRENTS (AMPS)) 

WRITE ( I PLACE * 30 ) 

30 FORMAT (3Xf 'LIFT ' r IX f ' PITCH ' f IXf ' AXIAL ' f 4Xf 2HFU f4Xf 2HFLf4Xf2HFPf 

+ 4X r 2HFS 1 4 X f 2HFA f 4X f 2HAU f 4X f 2HAL f 4X f 2HAP f4Xf2HASf4Xf 2HAA ) 

DO 40 1=1 f NUMBER 
J=<I-1>*16+1 

H=IDATA( J+4 )+IDATA(J+l)+IDATA(J+2)+IDATA(J+3) 

P= ( IDATA( J+4 ) +IDATA ( J-f 3 ) )-(IDATA(J+2)+IDATA(J+l ) ) 

A=IDATA( J+5 ) 

FS=FLOAT ( ( " 7777 ♦ AND ♦ IDATA(J+6))~2047)*0*01221 
AS=FLOAT ( < " 7777 « AND* IDATA(J+7) ) -2047 ) *0 ♦ 01221 
FP=FLOAT ( ( " 7777 ♦ AND* IDATA( J+8) ) -2047 ) #0 . 01221 
AP=FLOAT(< *7777. AND* IDATAC J+9) ) -2047) *0*01221 
FU=FLOAT ( ( *7777* AND* IDATA< J+10) ) -2047 ) *0 * 01221 
AL=FLOAT ( ( " 7777 * AND ♦ I DATA ( J+ 1 1 ) ) -2047 ) *0 . 01 221 
FL=FLOAT ( ( -7777*AND*IDATA< J+12) ) -2047 ) *0 * 01221 
AU=FLOAT( ( "7777* AND* IDATA< J + 13) ) -2047) #0*01221 
FA=FLOAT< < *7777* AND* IDATAC J+l 4) ) -2047 >*0*01221 
AA=FLOAT ( ( • 7777 * AND ♦ I DATA ( J+ 1 5 ) ) -2047 ) *0 * 01221 
C 
C 

ARRAY(If1)=H 
ARRAY ( I f 2 ) =P 
ARRAY ( I f 3 > =A 
ARRAY ( I f 4 ) =FU 
ARRAY < I f 5 ) =FL 
ARRAY(If6)=FP 
ARRAY < I f 7 ) =FS 
ARRAY ( If8)=FA 
ARRAY ( I f 9 ) =AU 
ARRAY ( I f 1 0 ) =AL 
ARRAY ( If 11 )=AP 
ARRAY ( I f 1 2 ) =AS 
ARRAY ( I f 1 3 ) =AA 
ARRAY ( I f 1 4 ) =AA+FA 
XASUM=XASUM+AA 
XFSUM=XFSUM*f FA 
C 

IF < I • EG ♦ 1 ) SIGSUMAX=0 ♦ 0 
SIGSUMAX=SIGSUMAX f AA f FA 
C 

TIME=FLOAT ( IDATA ( J ) ) 

WRITE(IPLACEf50) (HfPfAfFUfFLfFPfFSfFAf AUfALfAPfASfAA) 

50 FORMAT ( 316 f 2X f 10F6 ♦ 2 ) 

FASUM=FASUM+FA 
AASUM=AASUM+AA 
AXSUM=AXSUM+A 
40 CONTINUE 

ICOUNT = ICOUNT + 1 

AUERSC ICOUNT, 1 )=ICOUNT 

AUERS(IC0UNTf2)=XASUM/NR 

AUERS ( I COUNT f 3 ) =XF3UM/NR 

AUERS ( ICOUNT f 4 ) = < XFSUM+XASUM ) /NR 
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c 

c *************************************************** 

C ROUTINE TO CHECK FOR REJECTION 

C 

C 

C CALC AVERAGE OF THE GET OF "NUMBER" READINGS 

AVSUMAX^S I G5UMAX/ NUMBER 

C 

C CALC SUM OF SQUARED DIFFERNCES FROM AVERAGE 

SUMDIFFSQ=0 ♦ 0 
DO 700 IZ=1 f NUMBER 

SUMDIFFSQ=SUMDIFFSQ+ ( ARRAY ( IZ f 14 > -AVSUMAX ) **2 
700 CONTINUE 

C 

C CALC STANDARD DEVIATION OF SAMPLE 

STDEV1= ( SUMDIFFSQ/NUMBER ) ** . 5 

WRITE( IPLACEf#) ' St dev= ' f STDEV1 f ' average sum of 
+ axial currents=' t AVSUMAX 

c 

C CHECK FOR DATA POINT OUT OF RANGE BY USING 

C CHAUVENET'S CRITERION 

C THIS IS SET FOR A BLOCK OF TEN VALUES* 

C IF ' NUMBER ' IS CHANGED THEN CHANGE THE '1*96' 

DO 710 IZ=1fNUMBER 
IFLAG2 ( IZ)=0 

IF ( ARRAY < IZ f 14 ) ♦ GE • ( AVSUMAX* 1 * 96*STDEV1 ) ) IFLAG2 ( IZ ) = 1 
IF (ARRAY < IZf 14 ) ♦ LE* ( AVSUMAX-1 * 96*STDEV1 ) ) IFLAG2( IZ)-1 
710 CONTINUE 

c 

C RESET TO ZERO 

DO 715 IZ=1 f 1 4 
SUM< IZ )=0*0 
715 CONTINUE 

c 

C WORK OUT NEW AVERAGES 

ICNT=0 

DO 720 IZ=1fNUMBER 

IF ( IFLAG2< IZ ) * EQ ♦ 1 ) GOTO 730 

ICNT=ICNT+1 

DO 740 IZ2=1f13 

SUM < I Z2 ) =SUM ( IZ2 ) + ARRAY ( I Z r I Z2 ) 

740 CONTINUE 

SUM < 1 4 ) =SUM (14) + ARR A Y ( I Z f 8 ) + ARRAY ( IZ f 1 3 ) 

730 CONTINUE 

720 CONTINUE 

DO 750 IZ=1 f 1 4 
AV< IZ)=SUM(IZ)/ICNT 
750 CONTINUE 

c 

C STORE AVERAGES OF FQR'D AX CURf AFT AX CUR 

C SUM OF AX CURRENTS AND AXIAL POSITION 

AVERS ( I COUNT f 5 ) =AV ( 8 ) 

AVERS ( I COUNT f 6 ) =AV (13) 

AVERS < ICOUNT f 7 ) = AV < 1 4 ) 

AVERS ( ICOUNT f 8 ) = AV ( 3 ) 

c 

C PRINT OUT 'GOOD ' DATA AVERAGES 

WRITE(5f*) / Averages of "good" data*** ' 
WRITE(IPLACEf371) (AV( JJ) fJJ=1f13) 

371 FORMAT ( ' ' f3F6.0f IXf 10F6*2) 

WRITE (5 f 370) AV ( 14 ) 

370 FORMAT(' Average sum of axial currents 2 7 fF5* 2f ' amps ') 

c 
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C PRINT OUT FLAGS TO SHOW WHICH DATA WAS G00D = 0 f BAD-1 

WR I TE ( IPL ACE * 380 ) ( I FLAG2 ( J J ) r J 1 r NUMBER ) 

380 FORMAT ( ' REJ' f 10 < IX f 21 ) ) 

WRITE ( 5 1 * ) ' ' 

C 

C HAVE THE PROPER NUMBER OF SETS BEEN ANALYSED 

C IF NOT GOTO 953 

IF ( I COUNT ♦ NE ♦ ISETS ) GOTO 953 

C 

C PRINT OUT WHICH SET AND AVERAGES 

DO 189 ICT=1 r ICOUNT 

WRITE(IPLACEf360) < AVERS (ICTf JJ) f JJ=1f8> 

360 FORMAT < ' ' f 7 < F7 * 3 f 2X ) f F6 . 0 ) 

189 CONTINUE 

C 

DO 191 ICT=2f8 
AVBACH< ICT ) =0 ♦ 0 
DO 192 IZ=1 f ISETS 

AVBACH(ICT)=(AVBACH(ICT)+AVERS(IZfICT)/ISETS) 

192 CONTINUE 

191 CONTINUE 

WRITE ( IPLACE f 360 ) 9999 ♦ #1000 « f ( AVBACH< JJ ) f JJ=2f 8) 

c 

C PASS VALUES TO NEXT SECTION OF THIS PROGRAM 

AXAV=INT ( AVBACH < 8 > + * 5 > 

FAAV=AVBACH ( 5 ) 

AAAV = AVBACH < 6 ) 

c 

C RESET COUNTER TO ZERO FOR NEXT BATCH OF SAMPLES 

I COUNT = 0 

C END OF REJECTION ROUTINE 

C 

C INPUT ROUTINE FOR DATE AND FLAGS 

C 

WRITE ( 5 f 823 > 

WRITE ( 5 f 824 ) 

WRITE ( 5 f 825 ) 

WRITE<5f826) 

WR I TE < 5 f 827 ) 

WRITE ( 5 f 828 ) 

WRITE( 5 f 829 ) 

IGRP=IGRPF1 

WR I TE ( 5 f 834 ) I GRP 

834 FORMAT ( 7 THIS IS SET NUMBER' f3I ) 

823 FORMAT ( ' INPUT DATE EG 111286 ' ) 

824 FORMAT ( ' IF FIRST TEST ft WIND IS OFF THEN DATE=*START DATE" ') 

825 FORMAT ( ' IF LAST TEST ft WIND IS OFF THEN DATE='STOP DATE* ') 

826 FORMAT ( ' EG START 1112 OR STOP 1112 MAX 12 CHARS' 

+ f/f' THESE ARE USED TO FIND THE ZERO WIND CURRENTS ') 

827 FORMAT ( ' 'f/f' IF NOT START OR STOP AND ') 

828 FORMAT ( ' ' f!8Xf 'STING OUT THEN DATE=’NO DATE 1 EG NO 111286 ') 

829 FORMAT ( ' 'f18Xf 'STING IN THEN DATE= 1 YS DATE 1 EG YS 111286 ') 

READ< 5f 833 ) ( DATES (II) f 11 = 1 f 7) 

833 FORMAT ( 7A2 ) 

C 

C INPUT ROUTINE FOR BASE PRESSURES 

C RESET VALUES IN CASE START OR STOP RUN 

P1=0.0 
P2=0 • 0 
P3=0 • 0 

C ALSO RESET DYNAMIC PRESSURE IN CASE OF START OR STOP 

PA=0 ♦ 0 

C NO NEED IF NO STING HENCE GOTO 905 

IF < DATES < 1 ) #EQ ♦ ' NO ' ) GOTO 905 
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c 

800 

801 

905 

C 

820 

933 

910 

920 

840 

C 

C 

860 

870 

930 

C 

891 

C 

C 

C 

964 

C 


OR IF START OR STOP HENCE GOTO 933 
IF < DATES < 1 ) .EG • ' ST ' ) GOTO 933 
WRITE ( 5r 800 ) 

WRITE (5r 80 1 ) 

FORMAT ( 7 INPUT CENTRE MIDDLE & OUTER PRESSURES IN KPASCAL ') 
FORMAT ( / RELATIVE TO REFERENCE STATIC PRESSURE. ' ) 
READ<5r*)PlrP2rP3 
CONTINUE 

WRITE < 5? 820 ) 

FORMAT ( 7 INPUT REFERENCE DYNAMIC PRESSURE IN KPASCAL ') 

READ(5r*)PA 

CONTINUE 

IF(DATES< 1 ) .NE. 'ST' ,0R. C DATES ( 2 ). NE . 'OP' ♦ AND.DATES(2) .NE. 'AR' ) ) 
+GOTO 910 

WRITE ( 3 r 840) AX A V r FAAV r AAAV r PA r RTC r RP r PI f P2 r P3 r < DATES < 1 1 ) * 

+ 1 1 = 1 r 7 ) 

GOTO 920 
CONTINUE 

WR I TE ( 2 » 840 ) AXA V r F AAV r AAA V r PA f RT C r RP r P 1 r P2 ? P3 r (DATES(II) r 
+II=lr7) 

CONTINUE 

FORMAT < I8r2Xr 2 <F6.3r2X) rF7.3r 2XrF5*2r4(2Xr F7.3) r2Xr7A2) 


RUN AGAIN 
WRITE <5* 860 ) 

FORMAT ( ' DO YOU WISH TO RUN AGAIN (Y/N) ') 

READ (5*870) KEY 
FORMAT (Al) 

IF (KEY *EQ. ' Y' )GOTO 900 

WRITE(2r840)0r0.0r0.0r0.0r0.0r0.0r0.0r0.0r0.0r ' ♦ ♦ ' r 'EN' 

+ r'D 'r'FI'r'LE'r '..'r'**' 

REWIND 2 
CONTINUE 

READ < 2 r 840 ) AXAV r FAAV r AAAV r PA r RTC r RP r PI r P2 r P3 r 
+ < DATES (II) rll“l,7) 

WRITE(3r840)AXAVr FAAVr AAAVrPArRTCrRPr PI rP2rP3r 
+<DATES(II) r II=lr7) 

IF(DATES<1) .NE. ' . . ' . AND . DATES ( 7 ) . NE . ' . . ' )GOTO 930 
APPEND TO CURRENTS FILE? 

WRITE ( 5 r 891 ) 

FORMAT ( ' WRITE THE NEW DATA TO THE MASTER DATA FILE (Y/N) ') 
READ ( 5 r 870 >KEY 
IF(KEY. NE. 'Y' )GOTO 970 

IF YES THEN APPEND TO CURRENTS MASTER FILE 

CLOSE 4PR0G1.TEM FILE SO THAT CHANNEL CAN BE USED AGAIN 

CLOSE ( UNIT=2 ) 

REWIND 3 

OPEN CURRENTS MASTER FILE 

OPEN ( UNI T=2 r NAME= ' 4CURR . TEM ' r TYPE = ' SCRATCH ' r FORM= 'FORMATTED ' ) 
OPEN < UNIT=4 r NAME= ' 4CURR • DAT ' r TYPE = ' OLD ' r FORM= 'FORMATTED ' > 
CONTINUE 

READ FROM CURRENTS FILE & STORE IN TEM FILE 
READ( 4 r 840) AXAV r FAAVr AAAVrPArRTCrRPr PI r P2r P3r (DATES< II) r 
+II=lr7) 

IF ( DATES ( 1)«EQ •'.«') GOTO 965 

WRITE < 2r 840) AXAV r FAAVr AAAV r PA r RTC rRPr PI r P2 r P3 r ( DATES < II) r 
+ 1 1 = 1 r 7 ) 

GOTO 964 
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965 

C 


777 

C 

c 

c 


967 


C 


893 

970 

892 

890 


CONTINUE 

READ FROM 4PR0G1 ♦ STO f WRITE TO 4CURR*TEM INCLUDE END FILE 
READ ( 3 f 840 ) AXAV ? FAAV t A AAV f PA f RTC f RP f PI f P2 f P3 f ( DATES < 1 1 ) f 
FI 1 = 1 f 7) 

WR I TE ( 2 f 840 ) AXAV f FAAV f AAAV f PA f RTC f RP f P 1 f P2 f P3 f ( DATES ( 1 1 ) f 
FI 1=1 f7) 

IF ( DATES ( 1 ) ♦ NE. 7 ♦ ♦ 7 ) GOTO 965 
GOTO 777 

WRITE(3f840)0f0*0f0»0fC*0f0*0f0*0f0*0f0*0f0«0f 7 . * ' r'EN* # 

Ff 7 D 7 f 'FI 7 f 7 LE 7 f 7 ♦ ♦ 7 f 7 ♦ ♦ 7 
CONTINUE 
CLOSE ( UNIT=4 ) 

OPEN MASTER FILE AS NEW FILE THUS WIPING OUT OLD VERSION 
& ALLOWING NEW VERSION TO BE READ IN FROM TEMP FILE 
OPEN ( UNIT=4 f NAME= 7 4CURR* DAT 7 fTYPE= 7 NEW 7 f FORM 55 7 FORMATTED 7 ) 
REWIND 2 
CONTINUE 

READ ( 2 f 840 ) AXAV f FAAV f AAAV f PA f RTC f RP f PI f P2 f P3 f ( DATES ( 1 1 ) f 
FI 1 = 1 f 7 ) 

WRITE ( 4 f 840 ) AXAV r FAAV t AAAV f PA f RTC f RPf PI ? P2 f P3 f ( DATES ( 1 1 > f 
+ 1 1 = 1 f 7 ) 

IF ( DATES ( 1 ) .NE* 7 ♦ ♦ 7 )GOTO 967 

CLOSE DATA FILES WIPING OUT SCRATCH FILES 

CLOSE ( UNIT =4 ) 

WRITE(5f893) 

FORMAT ( 7 DATA HAS BEEN ADDED TO MASTER DATA FILE 4CURR.DAT 7 ) 

CONTINUE 

CLOSE ( UNI T=1 ) 

CLOSE ( UNIT=2 ) 

CLOSE ( UN IT =3 ) 

WRITE(5f892) 

FORMAT ( 7 THE DATA HAS BEEN STORED IN FILE 4PR0G1*ST0 7 > 

WRITE ( 5 f 890 ) 

FORMAT ( 7 THE LATTER FILE IS USED IN THE PROGRAM 4PR0G2 7 ) 
END 
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APPENDIX IV 


4PROG2 . FOR 


C 

C PROGRAM 4PR0G2 ♦ FOR 

C A #U# NEWCOMB 

C 

INTEGER AXAV 
DIMENSION DATES< 7) 

REAL LENGTH f MACH fNU 
C 

FSTART-O ♦ 0 
FST0P=0 *0 
ASTART =0 ♦ 0 
AST0P=0 * 0 

880 FORMAT < 7 '>4A10> 

OPEN < UNI T= 1 f NAME= ' 4PR0G 1 ♦ STO ' f TYPE* ' OLD ' f FORM* ' FORMATTED ' ) 

C ********** ALL PRESSURES IN KPASCAL ********** 

C ******** RELATIVE TO REFERENCE STATIC ******** 

OPEN ( UNIT=2 f NAME* ' 4PR0G2* STO ' t TYPE* ' NEW ' f FORM* 'FORMATTED' ) 

900 CONTINUE 

JJ*JJ+1 

WRITE( 5f 890) JJ 

890 FORMAT ( ' READING IN LINE# ♦ ♦ ♦ ' f 12) 

C READ IN DATA FROM 4PR0G1#ST0 

READ( 1 f 800) AXAVfFAAVf AAAVf PAfRTCf RPf PI fP2fP3f 
+ ( DATES (II)fII=1f7) 

800 FORMAT ( 18 f 2X f 2 < F6 ♦ 3 f 2X) fF7*3f2XfF5#2f4(2XfF7*3) f2Xf 7A2) 
IF(DATESd) .EQ* ' * * ' ♦ AND . DATES < 7 ) ♦ EQ* ' ♦ ♦ ' )GOTO 910 
C 

IF ( DATES ( 1 ) *NE« 'ST' #0R. ( DATES ( 2 ) ♦ NE# 'AR' ♦ AND ♦ DATES ( 2) . NE. 'OP' )) 
4-GOTO 945 

WRITE(5f880) 'CURRENT 'f 'OFFSETS 'f 'RESET' 

IF < DATES<2> *EQ ♦ ' OP ' ) GOTO 946 
FSTART=FAAV 
ASTART=AAAV 
GOTO 944 
946 CONTINUE 

FSTOP=FAAV 

ASTOP~AAAV 

944 CONTINUE 

C NO NEED TO CALCULATE VALUES AS CURRENT RESET RUN THROUGH 

GOTO 900 

945 CONTINUE 

C CONVERT CURRENTS TO FORCES 

CURROFFSET = ♦ 5* < FSTART + AS TART +FSTOP+ASTOP ) 

FGRCEFACTOR-O ♦ 04385 

AXFORCE=FORCEFACTOR* ( FAAV+AAAV “CURROFFSET ) 
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c 

IF ( PA ♦ NE ♦ PB ) GOTO 940 
WRITE < 5 f 804 ) 

WRITE(5f805) 

804 FORMAT ( ' DYNAMIC PRESSURES. 0 f RUN WILL CONTINUE ') 

805 FORMAT ( ' DATA WILL BE STORED AS OVERFLOWS ') 

C MAKE NUMBERS TOO LARGE FOR FORMATTED OUTPUT 

C OR SET TO 0,0 

CD=999 • 0 
REY = 0 ♦ 0 
V=0.0 
MACH=0 • 0 
PDYNTRUE=0. 0 
GOTO 960 


940 CONTINUE 

C 

C * ************************ 

C CALC TRUE TEST SECTION DYNAMIC PRESSURE 

DYNFACT0R=1 .017 
PD YNTRUE=PA*DYNF ACTOR 
IF < PA .GE.O.O)GOTO 930 
WRITE < 5 f 850 ) 

850 FORMAT < 7 DYNAMIC PRESSURE < 0.0 WILL TAKE MODULUS ') 

PDYNTRUE=PDYNTRUE*-1 
930 CONTINUE 

C 

C ***************** ******** 

C CALC ROOM DENSITY 

RHO=RP/ < 287* ( RTC+273 .15))* 1000 
C 

C * * * * * * * * ^ ^ * * * * * * * * * * * * * * ^ 

C CALC NU 

RECMU0=56085 . 3 

NU=( ( (RTC + 273.15)/288. 16) ** * 8 ) /RHO/RECMUO 
C 

C CALC BASE FORCE 

C IS STING IN PLACE? 

C IF NOT GOTO 905 

BASEF0RCE=0.0 

IF( DATES (1) *NE. ' YS ' )G0T0 905 
C STING IS IN HENCE PROCEED 

C ASSUMING QUADRATIC DISTRIBUTION 

C OF PRESSURES OVER THE BASE AREA* 

C ♦♦♦..HAVE RADII IN MILLIMETRES 

C TO MAKE NUMBERS BETTER 
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R1=0.0 
R2=l .9 
R3=3.8 
R4=4.0 
PI=3. 141593 


C IF RELATIVE TO REF STATIC 

C CP=P/PDYNTRUE+1-PDYN/PDYNTRUE 

C USE CP THEN CHANGE TO P LATER 


Pl=l .O+Pl/PDYNTRUE-1 .O/DYNFACTOR 
P2=l . 0+P2/PDYNTRUE- 1 . O/DYNFACTOR 
P3=l .0+P3/PDYNTRUE-1 .O/DYNFACTOR 
A=P1 

B=P2*R3*R3-P3*R2*R2-P 1 *R3*R3+P 1 *R2*R2 
B=B/ < R2*R3*R3-R2*R2*R3 ) 

C= ( P3-P1 ) * ( R3-R2 ) *R2*R3+P1 *R3* ( R3*R3-R2*R2 ) -R3* ( P2*R3*R3-P3*R2*R2 ) 
C=C/ ( R2* (R3-R2 ) *R3*R3*R3 ) 


C BASEFORCE ON CP = 

BASEFORCE=PI* < 3 . 0*C*R4#*4+4 . 0*B>ttR4*!)t3+6 . 0*A*R4*R4 ) /6 . 0/ 1 OOO . 0 

C BASEFORCE ON PRESSURES = 

BASEFORCE=BASEFORCE*PDYNTRUE 
905 CONTINUE 

C 

c a######*###***#*####*#### 

C CALC RESULTANT AXIAL FORCE 

C IF BASEFORCE POSITIVE AS THRUST 

RAXF=AXFORCE+BASEFORCE 

C 

C ************************* 

C CALC DRAG COEFFICIENT 

AREA=. 00047 

CD=RAXF/ ( ( AREA*1 000 ) *PDYNTRUE ) 

C 

C ************************* 

C CALC REYNOLDS NUMBER 

V= < PDYNTRUE* 1000*2 . 0/RHO ) ** . 5 
LENGTH=0. 184152 
REY=U*LENGTH/NU 

MACH=V/SGRT ( 1 . 4*287 . 0* ( 273 . 1 5+RTC ) ) 

C 

960 CONTINUE 

C 

C NOU HAVE ALL VALUES NEEDED 

C WRITE TO NEW DATA FILE 

WRITE <2, 810>CD» REY,V, MACH, PDYNTRUE, AXFORCE. BASEFORCE, 

+( DATES! II), 11=1, 7) 

810 FORMAT (F7 . 4 ,2X ,F7.0,2X,F5.2,2X, F6 .4,2X,F7.3,2(2X, F8 .5) ,2X, 7A2) 
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c 

c 

c 

c 

910 

820 


C 

C 

825 

826 
C 

C 


964 

C 


965 

C 


NOW LOOP BACK * READ IN NEXT SET OF DATA FROM 4PR0G1*ST0 
GOTO 900 

DATA READ IN FINISHED 
WRITE< 5f820) 

FORMAT ( 7 DATA FILE EMPTY NOW ') 

WRITE(2f810)O.OfO*OfO*OfO*OfO*OfO,OfO,Of 7 ♦ ♦ 7 f 7 EN' f 
+ 7 D ' f'FI' f 'LE' f 7 . . 7 f 7 ♦ ♦ 7 

APPEND TO DRAG FILE? 

WRITE ( 5 f 825 ) 

FORMAT < 7 WRITE THE DATA TO THE MASTER DRAG FILE <Y/N> 7 ) 

READ (5f 826) KEY 
FORMAT ( A1 ) 

IF<KEY*NE* 7 Y 7 )GOTO 970 
IF YES TO APPEND THEN 
REWIND 2 

OPEN DRAG DATA FILE 

OPEN ( UNIT=3 f NAME= 7 4CDRAG ♦ TEM 7 f T YPE= 7 SCRATCH 7 f FORM= 7 FORMATTED 7 ) 
OPEN ( UNI T=4 f NAME= 7 4CDRAG ♦ DAT 7 f T YPE= 7 OLD 7 f FORM= 7 FORMATTED 7 ) 
CONTINUE 

READ FROM DRAG FILE & STORE IN TEMP FILE UNTIL END FILE 
READ(4f810)CDfREYfUfMACHfPDYNTRUEfAXFORCEfBASEFORCEf 
+ ( DATES (II) fII=1f7) 

IF ( DATES (1)»ECK'»* 7 ) GOTO 965 

WRITE(3f810)CDfREYfVfMACHfPDYNTRUEfAXF0RCEfE<ASEF0RCEf 

+<DATES(II)fII=1f7) 

GOTO 964 
CONTINUE 

READ FROM 4PR0G2.ST0 ONTO TEMP FILE INCL* END FILE 
READ ( 2 f 810 ) CD f RE Y f U f MACH f PDYNTRUE f AXFORCE f BASEFORCE f 
+(DATES<II) fII=1f7) 

WRITE<3f810)CDfREYfVfMACHfPDYNTRUEfAXF0RCEfBASEF0RCEf 

+(DATES<II)fII=1f7) 

IF ( DATES ( 1 ) fNE* 7 « ♦ 7 )G0T0 965 

WRITE(3f810)0»OfO»OfO»OfO*OfO*OfO*OfO*Of 7 ♦ ♦ 7 f 7 EN 7 f 
+ 7 D 7 f 7 FI 7 f 7 LE 7 f 7 ♦ * 7 f 7 ♦ . 7 
CLOSE (UNI T= 4 ) 
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n o n 


OPEN DRAG FILE AS A NEW FILE. THIS WIPES OUT THE OLD DATA 
S ALLOWS OLD l NEW DATA TO BE READ IN FROM TEMP FILE 
OPEN < UNIT=4 > NAME= ' 4CDRAG . DAT ' » TYPE= ' NEW ' > FORM= ' FORMATTED ' ) 
REWIND 3 
967 CONTINUE 

READ ( 3 , 81 0 ) CD . RE Y » V > MACH » PDYNTRUE » AXFORCE . BASEFORCE » 

+ < DATES< II)»II = 1»7) 

WRITE < 4 » 810 ) CD » RE Y » V . MACH , PDYNTRUE > AXFORCE > BASEFORCE » 

+( DATES < II) » I 1=1 *7) 

IF ( DATES ( 1 > .NE. ' . . ' )GOTO 967 

C CLOSE DATA FILES WIPING OUT THE SCRATCH FILE ON 3 

CL0SE<UNIT=3) 

CLOSE ( UNI T=4) 

WRITE<5,839> 

839 FORMAT ( ' NEW DATA HAS BEEN ADDED ONTO DRAG FILE 4CDRAG.DAT '> 
970 CONTINUE 

CLOSE ( UNIT=1 ) 

CLOSE < UNIT=2 ) 

C 

WRITE<5.840) 

840 FORMAT (' NEW DATA HAS BEEN STORED SEPARATELY IN 4PR0G2.ST0 '> 
END 
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APPENDIX V 


4MSBS3.MAC 


+ 

.TITLE 

CONTR 

{ DECLARE PROGRAM TITLE 

t 

* 

« GLOBL 

CONTR t KEYED r CHARS t COMMAN 

.GLOBAL SYMBOLS 

f 

•MCALL 

TRANSL , ROTATE » F'HASEF , OUTPUT 

{AND REQUIRED MACROS 


.MCALL 

ADCOL f INTEG » START . TRANS2 . ADST 

• 

9 


.MCALL 

ADCOLS . TIME r HAST > ROLL 

9 

* 

.MCALL 

.EXIT 


m 9 GAIN 

BLOCK 



9 

hg: 

.WORD 

19000. t HEAVE 

GAIN 

pg: 

.WORD 

22000. {PITCH 


sg: 

.WORD 

19000. » SLIP 


ygs 

.WORD 

26000. { YAW 


ag: 

.WORD 

17000. $ AXIAL 


rg: 

.WORD 

4. .ROLL 



DEFINE HARDWARE ADDRESSES 


CLBS=170432 
0UT=167774 
ADSR=167770 
ADBR=ADSR+2 
{ CALL STARTUP ROUTINE 


{ CLOCK < B) STATUS 
» DIGITAL OUTPUT BUFFER 


CONTRJ START OSCI LL » CONMAN » CHARS » DATA {CALL START MACRO 

* 
f 

{SET-UP FOR TSX 
{ 


MOV 

♦LOKJOBfRO 

{LOCK JOB IN MEMORY 

EMT 

375 



MOV 

♦PRIfRO 

{SET-UP 

PRIORITY 

EMT 

375 



MOV 

♦SETSGfRO 

{SET-UP 

SINGLE CHARACTER MODE 

EMT 

375 



MOV 

♦SETNWfRO 

{SET-UP 

NO WAIT MODE 

EMT 

375 



MOV 

♦SETIfRO 



EMT 

375 




r 


INITIALISE CHARS STORE WITH VALUES DEFINED AT BOTTOM 


MOV <?#FACTOR » CHARS+26 . 

MOV 04STMAIN » CHARS+40 . 

MOV 04STREF » CHARS+50 • 

MOV *0. >CHARS+52. { 

MOV 04THRESH.CHARS+54. 


{ INITIALISE CHARS+26. 
* 

F 

9 

* 

9 
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f 


i*****************************#******$*t**$t****************t*t*****tttt 
f MAIN LOOP START » CALL KEYBOARD HANDLING ROUTINE 


start: 

JSR 

PC.KEYBD 

i CALLS KEYBD SUBROUTINE 


TST 

CHARS 

.ARE UE FINISHED ? (ESC) 


BEG 

SYNCHR 

t JUMP IF NOT 


RTS 

PC 

.RETURN TO FORTRAN IF YES 

i LOOP 
* 

SYNCHRONISATION 


9 

SYNCHRJ 

INC 

CHARS+138. 

i INCREMENT MASTER COUNTER 


BIC 

*176000 .CHARS+138. r REDUCE COUNTER TO 10 BITS 

clock: 

TSTB 

0#CLBS 

IIS CLOCK DONE COUNTING(TO OVERFLOW) 


BPL 

CLOCK 

f WAIT IF NO 


BICB 

#2OO.0#CLBS 

> CLEAR CLOCK DONE BIT 

2 

MOV 

#151777. 0*OUT 

1 SETS CHANNEL 14 D/A TO -2.5V APPROX. 

9 CHANGE RA TO 

DATA STACK 


r 

MOV 

SP. STACK 

.ARCHIVE SYSTEM STACK POINTER 

* 

MOV 

DATAfSP 

f SET UP DATA STACK 

9 

TIME 

$ STORE 

TIME <1> 

* 

9 

?*#* M***************************************************************** 

; RESET 

SOME VALUES? 



TST 

CHARS+12. 

} CHECK CHARS+12 C ? 3 

* 

BEQ 

AWN3 

» IF=0 BRANCH TO AWN3 


MOV 

CHARS+26. .FACTOR .RESET FACTOR C F 3 


MOV 

CHARS+40. .STMAIN PRESET STMAIN CM3 


MOV 

CHARS+50. .STREF 

.RESET STREF t R 3 


TST 

CHARS+52. 

» = 0.0? [ S 3 


BEG 

AWN4 

.IF - 0.0 GOTO AUN4 


MOV 

RF. STREF 

.RESET STREF WITH RF 


MOV 

#0. .CHARS+52. 

f RESET TO 0+0 

AWN4 : 

MOV 

CHARS+54.. THRESH .RESET THRESHOLD! T 3 


MOV 

*0.. CHARS+12. 

* RESET FLAG C ? 3 

AWN3 J f 
r 



BRANCHED TO HERE 

J EXTRA ROUTINE TO ALTER AX READING OF MAIN DRAG 

( SENSOR TO COUNTERACT DRIFT. 



MOV 

STREF *R0 

♦LOAD RO WITH INITIAL VALUE 

£ 

SUB 

RFrRO 

» SUB SENSOR FROM INITIAL 


CMP 

RO, THRESH 



BGT 

AWNI 

$DIFF>THRESH 


NEG 

THRESH 

. THRESH=-THRESH 


CMP 

RO. THRESH 



NEG 

THRESH 

.RESET THRESH TO +VE VALUE 


BLT 

AWNI 

fDIFFOTHRESH 

2 

JMP 

AWN2 

»-THRESH<= DIFF >=+THRESH 

9 

AWNi: 

AC1=X1 

AC2=X2 

LDCIF 

R0.AC1 

f AC1=R0=DIFF 


MULF 

♦FACTOR. AC1 

f AC1=FACT0R*DIFF 


DIVF 

♦TSD.AC1 

» AC1=FACT0R*DIFF/1000 


LDCIF 

AX.AC2 

» AC2=MAIN DRAG SENSOR VALUE 


SUBF 

AC1.AC2 

*AC2=MAIN-FACT0R*DIFF/1000 


STCFI 

AC2.AX 

.STORE AC2 AS AXIAL 

AWN2 J 9 

4 

f ~ 



JUMP HERE IF NO CHANGE NEEDED 
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I 




I 


{*********************************!M************************************ 
f STORE POSITION SENSOR INFO 


MOV 

FP»RO 




DAST 

RO 

f STORE 

FP 

<2> 

MOV 

FS»RO 




DAST 

RO 

f STORE 

FS 

A 

u 

V 

MOV 

AP»RO 




DAST 

RO 

f STORE 

AP 

<4> 

MOV 

AStRO 




DAST 

RO 

{STORE 

AS 

<5> 

MOV 

AXfRO 




DAST 

RO 

{STORE 

AX 

<6> 


****t********t***tt****************************t****t***************M 
COMPUTE MODEL POSITION AND ATTITUDE 


i AND 

COLLECT 

SOME CURRENT 

DATA 


? 

ADST 

8* 

{INITIATE FRONT ST/BD A/D 


A 

9 

MOV 

FP, R1 

{FETCH MAIN POSITION SENSOR 

SIGNALS 


MOV 

FS, R2 

{FP=FORWARD PORTfAS=AFT STARBOARD { ETC 


MOV 

AP, R3 

A 

9 



MOV 

AS fR4 

A 

9 


SHEAVE 

MOV 

FSIfRO 

fLOAD RO WITH FP 



ADD 

R2 , RO 

f ADD ALL FOR HEAVE 



ADD 

R3 , RO 

F 



ADD 

R4, RO 

F 



SUB 

♦27777 f RO 

f REMOVE HEAVE OFFSET (27777)******* 


MOV 

ROfHEAVE 

{STORE HEAVE POSITION 



ADCOLS 


{COLLECT FRONT ST/BD 

<7> 


ADST 

10. 

{INITIALISE AFT STARBOARD 


J PITCH 

MOV 

R4 , RO 

{LOAD RO WITH AS 



ADD 

R3,R0 

{ADD OTHER AFT 



MOV 

R1,R5 

{LOAD R5 WITH FP 



ADD 

R2 r R5 

{ADD FS TO FP 



SUB 

R5 1 RO 

{CALC PITCH 



MOV 

ROrPITCH 

{STORE PITCH ATTITUDE 



ADCOLS 


{COLLECT AFT STARBOARD 

<8> 


ADST 

11. 

{INITIALISE FRONT PORT 


fSLIP 

MOV 

RlfRO 

{LOAD RO WITH FP 



ADD 

R3,R0 

{ADD OTHER PORT 



SUB 

R2 r RO 

{SUBTRACT STARBOARDS 



SUB 

R4, RO 

* 

f 



MOV 

ROfSLIP 

{STORE SLIP POSITION 


? YAW 

MOV 

RIfRO 

{LOAD RO WITH FP 



SUB 

R2fR0 

{SUBTRACT OTHER FRONT 



ADD 

R4 9 RO 

{ADD AS 



SUB 

R3fR0 

{SUBTRACT OTHER AFT 



MOV 

ROf YAW 

{STORE YAW ATTITUDE 



ADCOLS 


{COLLECT FRONT PORT 

<9> 


ADST 

13. 

{INITIATE AFT PORT 


i AXIAL 

MOV 

AXfRO 

{LOAD RO WITH AX 



SUB 

STMAINfRO 

{REMOVE AXIAL OFFSET 



ASH 

♦3 f RO 

{AXIAL SENSITIVITY BY 4 



NEG 

RO 

{CORRECT SIGN 



MOV 

ROfAXIAL 

{STORE AXIAL POSITION 
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LOOP RATE INDICATOR FLIP 


MOV 


♦ 151777 , 0*OUT i SETS CHANNEL 14 D/A TO -2.5V APPROX. 


*******************************$****t********************************** 
COMPENSATOR WITH INTERLEAVED DATA AGUISITION 


INITIATE DATA AQUISITION 
ADCOLS 
ADST 15. 


f COLLECT AFT PORT <10> 

} FORWARD UPPER CURRENT A/D START 


f HEAVE 


PHASEF HEAVE, HT, 0.44, 8. 0,2. 5 
TRANSL CHARS, OSC ILL, 6 
INTEG CHARS/HINT , 14. ,HG 
MUL R3,R2 

MOV R2,H0 


t PHASE ADVANCE HEAVE POSITION 

;add in user demands 

t INTEGRATOR 

$ HEAVE GAIN. DROP LOW ORDER PRODUCT 
r STORE HEAVE DEMAND 


( COLLECT DATA AND INITIATE NEXT 



ADCOLS 

t COLLECT FU CURRENT 

<11> 


ADST 

12 ♦ 5AFT 

LOWER CURRENT 


t 

♦PITCH 

PHASEF 

PITCHfPT rO« 37 

♦PHASE ADVANCE PITCH 

ATTITUDE 


ROTATE 

CHARS ♦ OSC ILL ♦ 4 

♦ADD IN USER DEMANDS 



INTEG 

CHARS ♦ PINT >30* ♦ PG 

♦INTEGRATOR 



MUL 

R3rR2 

♦ PITCH GAIN 



MOV 

R2 9 PO 

♦STORE PITCH DEMAND 


♦ 

♦ COLLECT DATA AND INITIATE NEXT 




ADCOLS 

♦COLLECT AL CURRENT 

<12> 

£ 

ADST 

9. ♦ FORWARD LOWER CURRENT 


♦ 

9 SLIP 






PHASEF 

SLIPfSTf0*44f8*0f2*5 

♦PHASE ADVANCE SLIP POSITION 


TRANSL 

CHARS ♦ OSC ILL ♦ 4 

♦ADD IN USER DEMANDS 



INTEG 

CHARS ♦ SINT ♦ 36 • ♦ SG 

♦ INTEGRATOR 



MUL 

R3fR2 

♦SLIP GAIN 



MOV 

R2fS0 

♦STORE SLIP DEMAND 


r 

♦ COLLECT DATA AND INITIATE NEXT 




ADCOLS 

i COLLECT FL CURRENT 

<13> 


ADST 

7 • ♦ AFT 

UPPER CURRENT 


9 YAW 






PHASEF 

YAWfYT ♦0*37f6*0»3#1 

i PHASE ADVANCE YAW ATTITUDE 


ROTATE 

CHARSfOSCILLf6 

♦ADD IN USER DEMANDS 



INTEG 

CHARS ♦ YINT ♦ 48 ♦ ♦ YG 

♦INTEGRATOR 



MUL 

R3 ♦ R2 

♦YAW GAIN 



MOV 

R2f YO 

♦ STORE YAW DEMAND 


♦ 

» COLLECT DATA AND INITIATE NEXT 




ADCOLS 

i COLLECT AU CURRENT 

<14> 

4 

ADST 

14. > FORWARD AXIAL CURRENT 


f 

♦AXIAL 

PHASEF 

AXI AL ♦ AT ♦ 0 * 54 ♦ 12 ♦ 0 ♦ 1 ♦ 

6 ♦ PHASE ADVANCE AXIAL 

POSITION 


TRANS2 

CHARS ♦ 2 

♦ ADD IN USER DEMANDS 



INTEG 

CHARSfAINTfO* ♦ AG 

♦INTEGRATOR 



MUL 

R3fR2 

♦AXIAL GAIN 



MOV 

R2 ♦ AO 

♦ STORE AXIAL DEMAND 


9 

9 COLLECT 

DATA AND INITIATE NEXT 



* 

» 

ADCOLS 

i COLLECT FA CURRENT 

<15> 

9 

ADST 

6 ♦ fAFT 

AXIAL CURRENT 
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f MODEL OUT DETECTOR 


outon: cmp 

#6777f HEAVE 

f TEST HEAVE POSITION 

BPL 

MODIN 

? JUMP TO OUTPUTS IF MODEL 

* MODEL OUT 



CLR 

HO 

9 CLEAR ALL FINAL DEMANDS 

CLR 

PO 

fIF MODEL IS OUT 

CLR 

SO 

J 

CLR 

YO 

9 

CLR 

AO 

9 

CLR 

HINT 

9 INCLUDING HIGH ORDER 

CLR 

PINT 

? INTEGRATOR ACCUMULATORS 

CLR 

SINT 

F 

CLR 

YINT 

F 

CLR 

AINT 

F 


modin: 


f 

9 COLLECT DATA 

ADCOLS f COLLECT AA CURRENT <16> 

9 

i*********t****t*****t****tttt******t*****t***************tttt***tt****t 

f OUTPUT ROUTINES WITH INTERLEAVED POSITION SENSOR DATA COLLECTION 
9 INITIATE A/D 
f FRONT STAR/BD 


ADST 

F 

5f 

i START FS A/D 

F 

f FORWARD UPPER 

OUTPUT 


MOV 

HOfRO 

i FETCH HEAVE DEMAND 

SUB 

POfRO 

$ SUBTRACT PITCH DEMAND 

OUTPUT 

10000 

fCALL OUTPUT ROUTINE 

f FORWARD LOWER 

OUTPUT 


MOV 

POfRO 

» FETCH PITCH DEMAND 

SUB 

HOfRO 

* SUBTRACT HEAVE DEMAND 

OUTPUT 

* 

70000 

fCALL OUTPUT ROUTINE 

f 

ADCOL 

FS 

i COLLECT 

f AFT STAR/BD 

ADST 

3* 

.INITIATE AS A/D 

7 

fAFT UPPER OUTPUT 


CLR 

RO 

.SET RO TO ZERO 

SUB 

HOfRO 

.SUBTRACT HEAVE DEMAND 

SUB 

POfRO 

i SUBTRACT PITCH DEMAND 

OUTPUT 

110000 

fCALL OUTPUT ROUTINE 

fAFT LOWER OUTPUT 


MOV 

HOfRO 

.FETCH HEAVE DEMAND 

ADD 

POfRO 

J ADD PITCH DEMAND 

OUTPUT 

40000 


F 

ADCOL 

AS 

i COLLECT 

fFRONT PORT 

ADST 

4 « 

i INITIATE FP A/D 


f FORWARD PORT OUTPUT 
MOV SOfRO 

ADD YOfRO 

OUTPUT 50000 
f FORWARD STARBOARD OUTPUT 
CLR RO 

SUB SOrRO 

SUB YOfRO 

OUTPUT 100000 
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! . 

i 

ADCOL 

FP 

i COLLECT 

! ? AFT 

PORT 



i 

i ; 

ADST 

2. 

{ INITIATE AP A/D 

j f AFT 

PORT OUTPUT 


1 

MOV 

YO.RO 

» FETCH YAW DEMAND 

* I 

j 

SUB 

SOfRO 

{SUBTRACT SLIP DEMAND 

1 

OUTPUT 

30000 

{CALL OUTPUT ROUTINE 

j {AFT 

STARBOARD 

OUTPUT 



MOV 

SO'RO 

{FETCH SLIP DEMAND 


SUB 

YO.RO 

{SUBTRACT YAW DEMAND 

* 

OUTPUT 

60000 

{CALL OUTPUT ROUTINE 


ADCOL 

AP 

{COLLECT 

fAXIAL 



! 

j ; 

ADST 

0. 

{INITIATE AXIAL A/D 

I ' 

| f AFT 

AXIAL OUTPUT 


I 

MOV 

AO.RO 

{FETCH AXIAL DEMAND 


OUTPUT 

20000 

{CALL OUTPUT ROUTINE 


MOV 

AO.RO 

{FETCH AXIAL DEMAND 

4 

OUTPUT 

120000 

{CALL OUTPUT ROUTINE 

4 

ADCOL 

AX 

{COLLECT AXIAL A/D 


ADST 

1* 

{INITIATE REFERENCE SENSOR 

* 

9 



A/D CONVERSION 


MOV 

#4 .RO 


dbgb: 

DEC 

RO 

{ DELAY 


TST 

RO 



BPL 

DBGB 



f 

ADCOL RF t COLLECT SENSOR DATA 

9 

f DELAY FOR OUTPUT PORT TO PERFORM JOB BEFORE LOOP RATE FLOP 
JMP DBG7 

9 DELAY REMOVED 

dbg?: 

; LOOP RATE INDICATOR FLOP 
* 

9 

MOM #155777. <?#OUT {SETS CHANNEL 14 D/A TO +2.5V APPROX. 

4 

9 

i*************************************t********************************* 

i PRESERVE DATA ? 



TST 

CHARS+132. 

{IS < SET +VE (STORE) 

7 


BEQ 

NOSAMP 

{JUMP OUT IF NO 



DEC 

CHARS+132* 

{KEEP DATA AND REDUCE 

CYCLE COUNTER 


BR 

JUMPST 

{DONE HERE 


nosamp: 

SUB 

#32* fSP 

{DISCARD DATA 


T 

f RESTORE SYSTEM 

i STACKS f RETURN 

TO LOOP START 


r 

jumpst: 

MOV 

SPfDATA 

{ARCHIVE DATA STACK POINTER 


MOV 

STACKfSP 

{RESTORE SYSTEM STACK 


* 

• 

JMP 

START 

{BACK TO {START 



f 


i ************************* ********************************************** 
i DECLARE VARIABLES 


r 

9 

chars: 

♦ BLKW 

72. 

{STORAGE OF INPUT COMMAND DATA 


data: 

♦ WORD 

0 F 0 

{ADDRESSES OF DATA ARRAY 


stack: 

♦ WORD 

0 

{SYSTEM STACK POINTER 


oscill: 

♦ WORD 

0 

{ADDRESS OF OSCILLATOR ARRAY 
-67- 



{ 


heave; 

• WORD 

0 

slip; 

.WORD 

0 

axial; 

. WORD 

0 

pitch; 

• WORD 

0 

yaw; 

• 

• WORD 

0 

f 

fp; 

. WORD 

0 

fs; 

• WORD 

0 

ap; 

• WORD 

0 

as; 

• WORD 

0 

ax; 

• WORD 

0 

rf; 

* 

• WORD 

0 

t 

factor; 

• WORD 

0. 

thresh; 

• WORD 

0. 

* 

tsd; 

• WORD 

1000. 

stref; 

• WORD 

2849. 

stmain: 

f 

• WORD 

2630. 

f 

ho: 

• WORD 

0 

po; 

• WORD 

0 

so: 

• WORD 

0 

yo: 

• WORD 

0 

ao: 

41 

• WORD 

0 

9 

ht : 

• FLT2 

0* f 0 ♦ 

pt: 

♦ FLT2 

0. 9 0. 

st: 

• FLT2 

0. 9 0. 

yt: 

• FLT2 

0 ♦ 9 0 • 

at: 

* 

• FLT2 

0 • 9 0 • 

9 

hint: 

• WORD 

0 F 0 

pint: 

• WORD 

OfO 

sint: 

• WORD 

OfO 

yint: 

• WORD 

OfO 

aint: 

* 

f 

• WORD 

OfO 


TSX+ EMT BLOCKS 


lokjob: 

• BYTE 

13f 140 

pri: 

♦ BYTE 

0 9 150 

» 

• WORD 

127. 

9 

setsg: 

♦ BYTE 

Of 152 


♦ WORD 

'S 

setnw: 

• BYTE 

Of 152 


♦ WORD 

'U 

seti: 

♦ BYTE 

0 f 152 

• 

• WORD 

'1 

9 

• END 



{MODEL HEAVE POSITION 
» SLIP POSITION 
{AXIAL POSITION 
> PITCH ATTITUDE 
{YAW ATTITUDE 

{FORWARD PORT POSITION SENSOR 

{FORWARD STARBOARD 

{AFT PORT 

{AFT STARBOARD 

{AXIAL 

{AXIAL REFERENCE SENSOR 
• scalar 

{INITIAL THRESHOLD FOR 
{REFERENCE SENSOR 

{ TSD=iOOO» 

{INITIAL REFERENCE SENSOR VALUE 
{INITIAL MAIN AXIAL SENSOR VALUE 


{FINAL HEAVE DEMAND 
{PITCH DEMAND 
? SIDE DEMAND 
{YAW DEMAND 
{AXIAL DEMAND 

{INTERMEDIATE DATA STORAGE 
{FOR PHASE ADVANCERS 


INTEGRATOR ACCUMULATORS 
(HIGH AND LOW WORDS) 
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APPENDIX VI 


4MSBS.DAT 


THIS VERSION TO BE USED FOR DRAG TESTING WITHOUT STING 
OR TESTING WITH STING IF MODEL IS FORWARD OF THE STING. 


50 

i 

h 

make sure 

400 

10 

< 

take data 

87 

10 



87 

10 

< 


87 

10 

< 


87 

10 

{ 


0 

0 

0 

the end J 


**** 


THIS* VERSION TO BE USED IF STING IS INSIDE BASE CAVITY 
OR FOR PRESSURE DISTRIBUTION INVESTIGATION. 


50 

400 

87 

87 

87 

87 

0 


1 h 
10 < 
10 < 
10 -C 
10 -C 
10 < 
0 0 


make sure the integrators are on 
take data 


the end *#** 
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APPENDIX 3ZXL 


PHOTOGRAPHS 



THE BOLTZ BODY, SAMARIUM COBALT MAGNETIC CORE, THE PRESSURE MEASURING 
REAR BODY, THE DUMMY STING BODY 
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ORIGINAL PAGE IS 
OE POOR QUALITY 





isiii-ailpi 


\ 




■ : : j : W . : J 




CLOSE UP OF PRESSURE DISTRIBUTION MEASURING REAR BODY AND DUMMY STING 
REAR BODY 



TOP TO BOTTOM: THE STING SUPPORT SLEEVE, THE DUMMY STING, AND THE 


PRESSURE DISTRIBUTION MEASURING STING 
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THE AXIAL POSITION DETECTOR, NOTE THE LARGE WINDOW FOR THE MAIN DETECTOR 
AND SMALLER WINDOW FOR THE REFERENCE DETECTOR ORIGINAL PAGE IS 

"73- QB EOOR QUALITY. 
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